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The authors describe how acoustic phonons can directly tunnel through vacuum and, therefore, 
transmit energy and conduct heat between bodies that are separated by a vacuum gap. This effect is 
enabled by introducing a coupling mechanism, such as piezoelectricity, that strongly couples electric 
held and lattice deformation. The electric field leaks into the vacuum as an evanescent field, which 
leads to finite solid-vacuum-solid transmission probability. Due to strong resonances in the system 
some phonons can go through the vacuum gap with (or close to) unity transmission, which leads to 
significant thermal conductance and heat flux. 
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The heat flux of thermally excited photons from an 
ideal black body at a temperature T is given by the 
Stefan-Boltzmann law, which states that the flux is pro- 
portional to T 4 . The net flux between two black bodies, 
which are at different temperature and separated by a 
large vacuum gap d, is the difference of their individ- 
ual Stefan-Boltzmann fluxes. When a realistic emissiv- 
ity/absorption is considered this T 4 power law is altered, 
but the essential physics still remains the same. However, 
when d is smaller than the characteristic wave length At 
of thermal spectrum various near-field radiation effects 
start to play crucial role in the inter-body heat trans- 
port and new physics emerges (see Refs. [TJ [5] for a re- 
view). In this qxd < 1 limit (gy = 2tt/\t being the 
thermal wave vector) the heat flux is enhanced particu- 
larly by evanescent waves, as explained correctly first by 
Polder and Van Hove [3] . Recent advances in experimen- 
tal techniques have enabled near-field heat transfer mea- 
surements from (an down to 10 nm body distances. [4]-[6] 
In this Letter, we propose that at such distances a new 
type of evanescent field heat transfer mechanism due to 
acoustic phonons can exist. 

Even though acoustic phonons are the major heat car- 
riers in dielectrics, their effect on heat transfer through 
a vacuum gap has been considered to be negligible, be- 
cause they couple weakly to photons. Here we demon- 
strate by theoretical means that significant energy trans- 
mission and heat flux is possible if the acoustic phonons 
can induce an electric field, which then can leak into 
the vacuum [see Figs, [lja) and (b)]. Such a mecha- 
nism is provided, for example, by the density response of 
free carriers due to phonons or by the piezoelectric (PE) 
effect. Here we shall focus on the latter, which gives 
rise to a strong coupling between phonon induced ma- 
terial deformation and macroscopic electric fields. The 
solid-vacuum-solid transmission phenomenon described 
here can be thought of as an acoustic phonon tunneling 
through vacuum. 
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We consider a system in which two phonon black bod- 
ies 1/ and 3', which serve as fully thermalizing reservoirs, 
are connected to bodies 1 and 3, which serve as wave 
guides for the propagating acoustic phonons [Fig. [it a)]. 
The propagating modes hit the solid-vacuum interfaces 
and produce the evanescent electric fields of interest [Fig. 
[ljb)], which lead to non-zero solid- vacuum-solid energy 
transmission coefficient Xy = f° r incident mode 

7 = L,S. Mode indices L and S stand for longitu- 
dinal and transversal, respectively (the two transversal 
modes are not written explicitly). The energy transmis- 
sion probability 7^ 7 = 7^ 7 (g<i, 6) (from mode 7 to mode 
fi) is a function of the phonon polarizations, the abso- 
lute value of the incident phonon wavevector q and the 
angle of incidence 9. Now, following Ref. [7] the thermal 
boundary conductance G 7 arising from transmission of 
mode 7 can be defined as 
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where 73^ is the effective energy transmission coeffi- 
cient. The total thermal conductance and net heat flux 
are given by G — J2-y ^7 an d P = / T * GdT, respectively. 
Temperature dependency arises from phonon occupation 
number N(u} q ,T) = [exp{hu q /k B T) - l]" 1 and in T^ fl 
the temperature is buried into the thermal wave vector 
qx = 27r/ Ay = ksT/fiVj (At is the thermal phonon 
wave length and u 7 is the phonon velocity). For the 
sake of simplicity we have assumed a linear dispersion 



q in Eq. ( lb I . The parameter q c ~ 1/a is the 



Brillouin/Debye cut-off (a being the lattice constant). 
The bracket (• • •) stands for a solid angle average over 
half-space. Note that if 7^ = 1 and ^ — > 00, then 
= 1 and P is equal to the phonon black-body flux. 




FIG. 1: (Color online) Illustration of the phonon tunneling effect caused by evanescent electric fields, (a) Cross-sectional slice 
of area A illustrating a phonon incoming from a thermal path 1' at temperature T\ and hitting a solid- vacuum interface. The 
phonon carries an electric field, illustrated by +/- signs of polarization in between the wave fronts. The polarization induces 
an electric field into the vacuum gap 2 (some field lines illustrated). The field enables finite transmission over the gap into the 
thermal bath 3' (at temperature T3). The wave fronts of the reflected phonon are not shown, (b) A projection along the z-axis 
showing the spatial behaviour of the phonon waves (u) and of the evanescent electric field (E = Ee). The dashed curves depict 
the "reflected" evanescent field, (c) Illustration of the scattering of the incoming mode 7 into different propagating modes 
(L=longitudinal and S=transversal) . q 7 and p 7 are the propagation and polarization vectors, respectively, (d) The scattering 
matrix formulation of the problem (See Eq. [2b. 



The energy transmissions are calculated from the scat- 
tering matrices (S-matrix) Si (i = a, b) of the two solid- 
vacuum interfaces a and b. We define Si in such a way 
that it couples the amplitudes A^ v of propagating fields 
and exponentially decaying near-fields [see Fig. fljd)]. 
Here, j3 = ±, where + (— ) refers to left-to-right (right - 
to-left) propagating or decaying wave, a = 1,2, 3 is the 
material index (2 referring to the vacuum gap), and v is 
the mode/channel index (y — E refers to the evanescent 
channel). For example, for the interface a we write 
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The labels of the sub-matrices r a , t' al t a and r' a indicate 
the size of these matrices. We assume that in all materi- 
als there is only one channel that arises from evanescent 
electric fields. This fully covers the case in which we 
neglect all retardation effects and utilize the quasistatic 
approximation. Now, A® E is more conveniently related 
to the evanescent potential instead of the evanes- 
cent electric field Ee {Ee = — V$b)- The amplitudes 
A aE are coupled to the acoustic amplitudes if an acoustic 
phonon creates a periodic charge density or polarization, 
which then creates an evanescent field/potential due to 
a boundary. The retardation effects become important if 
the oscillation period of the surface polarization (arising 
from the acoustic phonons) is of the order of the time 
it takes light to make a round trip across the gap. This 
leads to cut-off energy £ c = hc/2d « 620 meVx(^m/d), 
below which our model is valid. Note that £ c is above 
acoustic phonon energies when d < 10 /im. 

For the full solid-vacuum-solid system we need to find 
the total S-matrix S — S a ® Sb that couples the ampli- 



tudes of the different solids. From S we specifically need 
the sub-matrix t, which describes the amplitude trans- 
mission. By solving S — S a <£> St, and taking into account 
the exponential factors arising from the finite distance d 
we find 
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where we have inverted the z-axis for Sb so that S a and 
Sb have a similar structure. Here the evanescent field is 
excited by the oscillating polarization perpendicular to 
the z-axis and, therefore, we have r\ = |sin#|. The energy 



transmission coefficients are given by 7^ 
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where a M7 is a factor that converts the amplitude trans- 
mission into energy transmission probability. We write 
the total energy transmission T~, = £\ T f ,~, in the form 
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where A^ = 

I = Im{ r^r^}. a M7 can be determined from the acoustic 
Poynting vector 
solid we have a u 



and in the simplest case of isotropic 
_ ^MT^gj^j^ 5 where p i is the mass 

1/2 

density of material i and {q^} z — (l — ( v fi/ v j) sm2 
is the z-direction propagation vector. 

Our next task is to calculate the energy transmission 
coefficient Ty in the case of two similar PE crystals sepa- 
rated by a vacuum gap. In such a piezoacoustic system, 
we are dealing with coupled acoustic and electromagnetic 
fields. Within the quasistatic model, the relevant field 
quantities (see Table [l]) and their couplings are defined 



TABLE I: Field quantities and constants. 
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Size 


Name 
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6x1 


Stress 




3x6 


Piezo tensor 
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lxl 


Potential 
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6x6 


Stiffness tensor 


u 


3x1 


Lattice displacement 


D 


3x1 


Electric displacement 
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3x3 


Dielectric constant 



by [5] 
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where V = [d/dx,d/dy,d/dz] T and V u is the 
displacement-to-strain operator (6x3 matrix). The 
form of the operator V u follows from the relation t a p = 
\(du a /dj3 + dup/da) (a, (3 = x,y,z), which couples the 
six symmetric strain components e a ,g to the components 
of displacement u. Note that Eq. (5c) is basically a Pois- 



son equation and the source term V eV ' u u gives rise to 
the phonon induced evanescent electric fields. The S- 
matrices Si are solved from a boundary condition equa- 
tion that is obtained by requiring the continuity of the 
potential $, the normal component of electric displace- 
ment D and the normal component of stress T at the 
interfaces. The boundary condition equation and the re- 
sulting Si are given in the supplementary material [S]- 
We calculate T 7 [Eq. Q] by solving the S-matrix nu- 
merically. We adopt material parameters that are close 
to that of ZnO [5] with the simplifying approximation 
= ^3j^i3e33, where e 33 = 1.3 C/m 2 . Furthermore, 
we assume that the acoustic properties are isotropic and 
that the PE stiffening can be neglected for the propa- 
gating modes. These assumptions have very little ef- 
fect on the angular averaged quantities of interest (T^^ 
and G 7 ). We use isotropic values en = 209.7 x 10 9 
N/m 2 ,c 44 = 42 x 10 9 N/m 2 and {e}^ = 10e <%. The 

phonon velocities are given by vl — \Jc-xxl p — 6119.3 
m/s and vg = J C44J ' p — 2738.6 m/s, where the mass 
density p = 5600 kg/m 3 . 

Figure [2] shows (qd, 9) for phonon tunneling from a 
PE material to another across a vacuum gap. There are 
two local minima in Ts'- one at 9 ss 25° and another one 
at 9 40° . The position of the latter minima can be re- 
lated to the phase change of the evanescent electric field 
in the vacuum-solid reflection. The former minima is 
the critical incidence angle where the quasi-longitudinal 
surface mode is exited. Most important effects for the 
energy transmission are the strong resonance features for 
both modes: for some (qd, #)-values T 7 is very large or 
even equal to unity. The resonances arise from the multi- 
ple reflections of the evanescent field in the vacuum gap, 
which leads to the [(e 27iqd - R) 2 + I 2 ] ~ l factor in T 7 [see 




FIG. 2: (Color online) Logarithmic contour plot of energy 
transmission coefficients 7^ (7 = L, S) between two PE bodies 
as a function of normalized wave vector qd and the angle of 
incidence 8. White regions have Ty < 10 -8 . The lower panels 
show log-log blow-up of the small-# large-gd region. 



Eq. d4||]. The elements r' a b are dictated by the vacuum- 
solid Doundary conditions and their real parts or moduli 
are not limited below unity. Thus we may have R > 1 
and / <IC 1, which leads to a sharp resonance peak at 
2rjqd — lni?. Indeed, for 7 = L, we have R ~ 1.5 for all 
incident angles and this is the origin of the sharp max- 
imum trajectory in Th{qd,9). For 7 = S, R > 1 only 
for angles 9 < 68°, and above this threshold no sharp 
resonances exist. For small 9 at qd 3> 1, 7~l = 1 at the 
resonance, whereas the amplitude of the resonance peak 
in Ts decays as a function of qd (see the lower panels of 
Fig. [2]). Note that the phonon resonant tunneling here 
has a striking similarity to the resonant photon tunneling 
described, for example, in Ref. [2]. 
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Inserting the calculated T-y into Eq. 
ing numerical integration, we obtain 7^ e 
results are presented in Fig. [3j We observe that at low 
temperatures when g^ti < 1 the effective transmission 
Tg^ is ~ 15% of the unity transmission resulting in rel- 
atively large Gs, which is of the order of the maximum 
possible thermal conductance allowed by Eq. (flaj). The 



large Tg^ at low temperatures follows from the broad- 
ened resonances and from the large transmission ampli- 
tude for the S'-mode for large 9 and small qd (see Fig. [2| . 
For 7 = L we have sharp resonances leading to a small 
effective transmission and to Gl <C Gs at low tempera- 
tures. 

Next, we want to investigate the dependency of T^'> 
on qxd and the contribution of different phonons in the 
(qd, 9) phase space to the heat transport. We divide 

into small angle 7^ e g<g th and large angle 7^g>e contri- 
butions: 
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FIG. 3: (Color online) Heat transport between two PE bodies: 
effective transmission T^' (the inset) as a function of qrd 
and the interface thermal conductance G 7 (the main part) as a 
function of temperature. The curves are obtained from the 7^ 
of Fig. |2| and Eq. |l| (we have set q c /<lT — > oo ). The dashed 
curves in the inset are small angle contributions g<e with 
6 t h = 3°. The dot-dash curve identicates (grd) -3 slope. For 
the y-axis units of the main Figure we have 27r 2 ti 7 fe_B/30d 3 = 



Y-, x (100 nm/d) 3 W/Km 2 , where Y L{S) = 55.6 (24.9). 



where the angular integrals of the different contributions 
are limited by the threshold 9 t h, which is chosen in such 
away that major contribution to 72 e<.Q th comes from the 
resonance trajectory at qd > 1. The exact choice of 9 t n 
is, therefore, somewhat arbitrary and in the numerical 
calculations we define 9 t h — 3°, whence 9 < 9 t h roughly 
correspond to the phase space of the lower panels of Fig. 
[2j At the high temperature limit (qrd ^> 1), small and 
large angle effective transmissions can be written as[H] 



r 



eff 



n+1 



y,v<Vth 



^ 15 [miZ(0)] 
^ ^ (2q T d) 2+n 

q T 



,1 



1th 



r 



eff 



15 1 



.; \ , dyy 2 Tj{y, 



(7a) 
,(7b) 



where F n (^,6 th ) 



9c /q T 



dx 
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with x t h 



\nR(9 th )/(2q T dsm9 th ), / M7 (0) = A^(9)/I(9) and 
fJL(0) is the first non-zero derivative of f lll (9). Due to 

the decay of the resonance for the S-mode, Ti 



■eff 

s,e<e th 



is 



very small. This means that practically all contribution 
to Tg and Gs at any temperature comes from the large 



■eff 
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At high temperatures, T^Ll 



7,8>e th 



is 



angle part Tl 

given by Eq. (7b), which explains the (qrd) J fall-off of 
T s e// and the saturation of G s G s s at in Fig. ^ The 



saturation value G s s at 



4TT 2 d 3 

follows from Eqs. (7b) and (Jlaj). For the L-mode the 
high temperature T^ n has also a large contribution aris- 
ing from 7^g<0 . This contribution eventually exceeds 
the large angle part and, as a result, there is no satu- 
ration in Gl- Close to 6 = 0, Tl 9<6 ls dominated by 
and b y setting ^ -> oo in Eq. 

exp(- 



q Q ° d dyy 2 Ty(y,e) 
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eff 



' L,e<e th ^ kh-i ^j /Uh=x TT 
{qrd) 5 ^ 2 behavior. It should be noted that the temper- 
ature dependency of the large angle contribution is not 
affected by the ratio q c /qr- As q^ 1 is of the order of 
the lattice constant a, Eq. (7b) is valid if d > a. The 
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we find 



resembling a 



small angle contribution depends on q c /qT- For example, 
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Thus, Gl also saturates when the cut-off q c is exceeded. 

At this point we summarize our findings: We have for- 
mulated the thermal boundary conductance of a solid- 
vacuum-solid system [Eq. Q] using an acoustic phonon 
energy transmission probability [Eq. Q] given by scat- 
tering matrix [Eq. (2l)], which couples the different solids 
by an evanescent channel produced by phonon induced 
electric fields (Fig. [lj. These fields can lead to a sig- 
nificant energy transmission probability for some parts 
of the wave vector-angle of incidence (qd, 9) phase space 
(Fig. |2]). The transmission exhibits resonances, which 
arise from the multiple reflections of the evanescent field 
in the vacuum gap. The L-mode shows a strong and 
sharp resonance trajectory across the (qd, #)-plane with 
a unity energy transmission probability. The resonance 
trajectory approaches 9 — as qd — > oo. The S-mode 
exhibits similar features. However, for a small qd and 
a large 9 the resonance is broadened and in the qd ^> 1 
limit the energy transmission decays rapidly. At low tem- 
peratures (qrd < 1), the broadened resonance at small 
qd leads to a large contribution from the S'-mode to the 
thermal conductance (~ 15% of the maximum possible 
thermal conductance defined by unity transmission). The 
small qd /large 9 part of the phase space dominates the 
S-mode heat transport leading to a saturation of the 
thermal conductance of the S-mode at high temperatures 
(Fig. [3J. If the thermal phonon wave length At is suf- 
ficiently smaller than the lattice constant, the L-mode 
shows no such saturation. This is due to the strong res- 
onance in the transmission that persists even at qd ^> 1 
and, therefore, allows short wave length phonons to go 
through the vacuum gap. Note that our method to cal- 
culate the thermal conductance is essentially similar to 
that adopted in (ballistic) thermal boundary resistance 
calculations [7], which means that the thermal excitation 
of local acusto-electric modes related to the solid- vacuum 
boundaries is not taken into account. Here, the bound- 
ary effects and the evanescent fields are solely due to 
radiation of bulk acoustic phonons. One can also explic- 
itly take into account thermal excitation and couplings 
of the local modes, expected to enhance the interbody 
coupling and the thermal conductance. This approach, 



5 



typically adopted in photon near- field studies [TJ will 
be left for future investigations. 

Acoustic phonon heat transfer through vacuum can be 
experimentally verified with direct thermal conductance 
measurements, using piezoelectric materials. The trans- 
mission of acoustic energy across a vacuum gap at long 
wave lengths can be investigated by measuring the trans- 
mission of a single acoustic wave by utilizing, for example, 
surface acoustic wave devices. Finally, we note that the 
effects described in this paper can also contribute to the 
thermal conductivity of polycrystalline piezoelectric ma- 
terials and percolation systems consisting of piezoelectric 
particles. Furthermore, as free electrons couple strongly 



to electric fields, a similar heat transport effect, that is 
discussed in this Letter, can also occur between metallic 
and piezoelectric bodies. 
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SUPPLEMENTARY MATERIAL 

i 



I. DERIVATION OF THE SCATTERING 
MATRIX IN A PIEZOACOUSTIC 
SOLID- VACUUM PROBLEM 

In the quasi-static approximation, the piezoelectric 
constitutive relations for stress T, electric displacement 
D and electric potential <£, generated by a displacement 
u, can be written as [5] 



where we have excluded the amplitude-phase factor 
Ae~i q r . The evanescent electric field is described 
through evanescent potential <&£ = A E e~ : > qp " r+q ' !Z (q p = 
Qx + Qy), f° r which similar field relations are given by 



T E = e T V<S> E = -jqe T L E 



-eV$ B = jqsL E , 



(14a) 
(14b) 



with 



T = e J V$- 
D = -eV$ 



cS/ u u 
- eV„it 



V T eV$ 



V T eV u it. 
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(8b) 
(8c) 



The electric field E is given by E = — V<&. Tensors e, 
e and c are dielectric permittivity, piezoelectric coupling 
and stiffness tensors, respectively. Above Equations are 
written following the abbreviated tensor subscript nota- 
tion [5j. The field quantities are are summarized in the 
Table I of the main article. Here, we have defined oper- 
ators V and V„ according to 

V T = [d/dx d/dy d/dz] (9a) 

d/dx d/dz d/dy 
d/dy d/dz d/dx . (9b) 
d/dz d/dy d/dx 



For displacement u with a known wave vector q, ampli- 
tude A and polarization vector p 



u = Ape jq r , 

the above operators can be written as 

V = -jqL 
V u = -jqL Ul 



(10) 



(11a) 
(lib) 



where L and L u are defined with the components l x , l y 
and l z of the unit propagation vector parallel to wave 
vector q: 



i x o o o i z l v 
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ly I 
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ly ij. 







Re-writing the constitutive relations for u 
with L and L u , we get 
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L T eL 
L T eL u p 
L T eL ' 
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eL u p 



(12a) 
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(13a) 
(13b) 
(13c) 



L% — [l x ly Ie] 



(15) 



where l e = q e /jq = ± \q p \ /jq = ±l/V2j. 

Boundary conditions at free piezoelectric boundary be- 
tween materials 1 and 2 (vacuum) require the absence of 
stresses and the continuation of <£> and the normal com- 
ponent of D. Thus, we have boundary conditions 





Don 



2«2 



= <J?9 



(16a) 
(16b) 
(16c) 



Txni 

where h a (a = 1,2 and |) 
tors of the boundary pointing into material a. Eq. ( 16a 
involves tensor form of T 



1) are the normal vec- 



Ti T 6 T 5 
U T 2 T 4 
T 5 T 4 T 3 



(17) 



Now, we include the notation of the main article, 
which writes the amplitudes as A^ avl where f3 = +(— ) 
denotes the left-to-right (right-to-left) propagating (or 
evanescent) wave, a denotes the material, and v de- 
notes the wave / field mode (L=longitudinal, S^shcar, 
_E=evanescent). From Eqs. (13), (14) and (16) we find 
the boundary condition matrix equation 



[T+n x ] [T-^] [T-^i] [T^ E n x ] 



D + rt 1 
iy l7 7l l 

17 
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A-XE. 



[0] 



^2E ll 2 
V 2E 
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'lE""i 
®2E 



A 2 E. 



(18) 



Here, A^ is the amplitude of an incident propagating 
mode 7 = L, S in material 1. The acoustic field related 
elements Tg„, and $g„(i/ = S,L) are given by Eqs. 

(13) and the evanescent field related elements Tt p and 



D aE are give by Eqs. (|14|). The elements <& p aE are de- 



fined by the continuity or the potential and, therefore, 
^ aE = 1. The elements in square brackets (for example 
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[T, B ni] ) are not scalars: here dim {[a;]} = 2x1. From 
Eq. (|18|) we obtain the scattering matrix equation 



Ail 
A\s 

A + 

A 2E 



= s 



A-iE. 



(19) 



where the scattering matrix S is defined by 



S 



[T-Si] [T^ni] [T^n,] [0] 



D 



is"i 



D 



lE'n 



D 



is 



4> 



IB 



2E' L 2 



-D+nT D 



'i 7 «i 
17 



[0] 

2E n 2 



2 A' 



(20) 



This equation is utilized in solving the S-matrix and the 
energy transmission in the main article. Note that all q- 
vectors v are proportional to the incident mode wave 
number q = \q^y\ and, therefore, S does not depend on 
q. Furthermore, as incident mode 7 is a parameter of S 
we must actually consider S as a sub-matrix of the total 
S-matrix. 



where #0 is defined by the equation qd = (2r] ) 1 In R(9 Q ) 
(i] = sin6>o)- Now for T^^ <e we can write 

Qc/q T 




I 2 cos 9 yR(8 ) 8<8 th 



7T 4 dqx 



Qo/q T 
x / dx- 



x 3 irA^^o) sin 6*o 



e x -l J(0 O ) 2 cos 6»o 



(25) 



Here x t h = ^R(9th)/(2qTdsm9 t h)- Let us further ap- 
proximate (note that sin6> = ln(R) /2xdqr) 



sin 6» _ lni?(6» ) 
cos 6*o 2xdqx ' 



(26) 



when we obtain 



II. 77 // AT g T d> 1 LIMIT 



In the main article we divide the effective transmission 
T* ff into small angle and large angle 7^f 9 f >8th 



contributions: 



'7 



// _ H-eff 1 T-e// 

- 'j,8<8 th 'i,8>8 th : 



(21) 



where the angular integrals of the different contributions 
are limited by the threshold 9 th- We will describe the 
derivation of the high temperature (qrd 3> 1) analytical 
formulas of T^tCa and 7T ,JLa below. 



A. Small angle contribution 7~^g{ C Q th 



t; 



E Tf 



15 1 



" "-" ' ^^{2dq T f 



' /qT dx Xirf^(6 )\nR(9 ) 
e x -l 



, (27) 



where / M7 (0o) = A^ 1 {9q)/I(9q). Next we perform a Tay- 
lor series expansion for f(6 ) = A tll (9 )/I(9 ) close to 
9 = 0: 



fnWo) = A^(e o )/I(0 o ) 



anF(0) 



/ ln(77(0)) \"/ M " 7 (0) 
\ 2dqxx J n! 



(28) 



The energy transmission coefficient is given by (see the 
main article) 



where / n (0) is the first derivative that is non-zero. We 
find 



T, 



^2-qqd 



(e 2r >qd _ ny + p 



Ea 



M7> 



(22) 



where = a M7 



{*iM*.} 7 . R = M^} and 7 = 



Im-jV^}. We use approximation (77 > 1, 7 <C 1) 



S{e 2riqd -R), 



(23) 



where S(x) is the Dirac delta function. This leads to 
expression 



7/J7 



„2 W ^ 5{9-6 Q ) 
I 2cos9 yR(9 ) 



(24) 



i,e<e th 



E 
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n+1 



„ ^ 3 (2 9T 7) 2+n 

/" 7 (0) p ( Qc „ , 
n\ q T 



where F, 



(-stew,) = f c/9r *;t . . 

,\q 1 1111 J x t h exp(x) — l 



(29) 



B. Large angle contribution 7~ e £{ 



At large angles only long wave length phonons partic- 
ipate, i.e., the phonon distribution can be approximated 
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by N(cu q ,T) = [cxp(fnu q /k B T) - l]" 1 w k B T/huj q and 
we find 



15 1 



q c d 



/ dyy%(y,6)} . (30) 



e>e 



th 



